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Abstract 

The electric and magnetic dipole moments of the rho meson are calculated 
using the propagators and vertices derived from the QCD Dyson-Schwinger 
equations. Results obtained from using the Bethe-Salpeter amplitude stud- 
ied by Chappell, Mitchell and Tandy, and Pichowsky and Lee, are compared. 
The rho meson EDM is generated through the inclusion of a quark electric 
dipole moment, which is left as a free variable. These results are compared to 
the perturbative results to obtain a measure of the effects of quark interac- 
tions and confinement. The two dipole moments are also calculated using the 

phenomenological MIT bag model to provide a further basis for comparison. 
13.40.Em, 14.40.Cs, ll.10.St, 12.38.Lg, 12.39.Ba 
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I. INTRODUCTION 



The electric and magnetic dipole moments of the rho meson, the lowest lying vector 
particle state are calculated using the semi-phenomenological n-point functions of Quantum 
Chromodynamics (QCD) obtained in the Dyson-Schwinger/Bethe-Salpeter framework. The 
rho meson is studied to provide some insight into the effect of QCD and confinement on the 
electromagnetic moments of hadrons, especially the question of the existence of a non-zero 
EDM for the neutron and the applications this has to the study of CP violation. In practice 
one is interested in these parameters for the nucleon and other hadrons, where they may be 
measured. But the study of the rho allows us to obtain insights without venturing into the 
complications of the three body system in QCD. 

The Dyson-Schwinger equations are a series of coupled integral equations relating the 
n-point functions of QCD to one another. They provide a continuum method of calculation 
in the strong coupling regime of QCD and as such are able to relate the Greens functions 
of QCD to hadronic properties. For a review of this topic see The use of the Dyson- 
Schwinger and Bethe-Salpeter equations to calculate hadronic properties has been applied 
to pion observables |§-[5[], kaon electromagnetic form factors |J, p-u mixing and the 
anomalous 771-* — ► mr form factor |J. 

This calculation makes use of a model form for the quark propagator which ensures con- 
finement by having no Lehmann representation, preventing free quark production thresholds, 
i.e., the quark cannot go "on mass-shell". A phenomenological form for the rho amplitude, 
which models the momentum distribution of the quarks inside, and the finite size of, the 
rho, is used, as well as the Ball-Chiu/Curtis- Pennington quark-photon vertex which ensures 
that both the Ward-Takahashi Identity and multiplicative renormalizability are preserved. 
The CP violating electric dipole term is introduced via a quark EDM which is left as a free 
parameter which can be varied. 

Results obtained from using the Bethe-Salpeter amplitude studied by Mitchell and Tandy 
0, Chappell ||, and Pichowsky and Lee RIQ| , are compared. Perturbation theory results 
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Tl|| , which ignore the dressing of the vertex functions and propagators, are also included to 
observe the effect of using a quark propagator that correctly incorporates confinement and 
of modelling the dressed quark core of the rho. 

A phenomenological model, the MIT bag model is also examined. This model incorpo- 
rates relativistic, point-like confined quarks but ignores one-gluon exchange and centre of 
mass motion. Both the magnetic and electric dipole moments are calculated in the bag model 
providing a basis for comparison between the Dyson-Schwinger/Bethe-Salpeter framework 
and a purely phenomenological model, one that has been well studied over the years. 

II. THE DIPOLE MOMENTS 
A. The pp7 Vertex 



The form for the CP conserving pp7 vertex is 12 



d:(p)T^ a d°(p') = d:( P )[( P + P '),(-g ua £(q 2 ) 

+quqaQ(q 2 )) + {g^Qa ~ 9^q v ) M(q 2 )}d° (?/), (1) 



where d a/ 3 are the rho polarisation vectors, 

k a kg 

d a p{k) = g af3 - (2) 
The form factors S, Q and M. reduce in the limit as q 2 — > to 
S(0) = 1, 

A4(0) = fi, the magnetic moment in units of e/2m 
2 



Q(0) = — - (Q + fi — 1) , where Q is the quadrupole moment in units of e/m 
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B. Projecting Out Charge and Dipole Terms 



The integrals for the rho dipole moments that will be derived later, equations (10) and 
(39), are too difficult to perform analytically. This means that to isolate relevant terms they 
will have to be projected out from the integral before it is calculated numerically Defining 



(3) 



and with the projection operator 



pr = ^( P )^(-/ ff ')Cb') 



(4) 



where P M = (p + the following result is obtained, 



M V al j,i3 



125(0). 



(5) 



To project out the magnetic dipole term we use the following projection operator, 



P^ = d\,{p) 



6p 2 



dUp) 



(6) 



which is such that 



occuP-ir lim<j 2 ^0 
% Van/3 — > 



4M(0). 



(7) 



C. Quark Contribution to the Magnetic Dipole Moment 

All calculations below, except those involving the bag model, are carried out in Euclidean 
space where = S^, 7m = 7^ an d {7^,7^} = 25^. 

The impulse approximation to the quark contribution to the magnetic dipole moment is 
given by (see Fig. |B|) 

M fi = £ e a (p)e^q)e^p')eQ &avon J^. (8) 

flavour=u,d 

Qflavour is the charge of the quark interacting with the photon. The u and d-quarks are 
treated as identical except for their charge and so the contribution becomes 



M ft = e a (p)e^q)e^p')el a ^, 



(9) 



W = (-l) / ^ I tr CFi ,[r^(fc;p-g)S(A;_ + )zr^ + ^(A; ++ ,A;_ + ) 

xS(k ++ )T^k + q/2;p)S(k + ^}. (10) 

Where fc a/3 = k + f + T p refers to the rho meson amplitude and Y^ C+CP is the dressed 
quark-photon vertex, and S(p) is the dressed quark propagator for a quark of momentum 
p, which will all be discussed in detail below. 



D. Quark Propagators 

The general form for the solution to the quark propagator Dyson- Schwinger equation 



is 



S(p) = -i i>ay{$ 2 ) + a s (p 2 ) 

= (i m(p 2 ) + B(p 2 )y l . (ii) 

A model form for the propagator is given by |3],f|] 

fl - e~ 2 ^ 2 A , (12) 
x + m z v ' 

, 2{x + m 2 ) - 1 + e -2(-+™ 2 ) 
avi?) = -± / _ mC m e 2x , 13 

where x = p 2 /2D, av(x) = (2D)av{p 2 ), &s = V^DcFsip 2 ) and rh = m/\/~2D, D is a 
mass scale. (A = 10~ 4 is chosen to decouple the small and large spacelike— p 2 behaviour in 
equation (12); i.e., to allow for 6 to govern the ultraviolet behaviour and 62 the infrared.) 
The parameters C m , in, bo, ■ ■ ■ , 63 are || 
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= 0.0 


Cm=0 


= 0.121 


m 


= 0.00897 




= 0.131 




= 2.90 


b 2 


= 0.603 


h 


= 0.185 



(14) 



with the mass scale D = 0.160 GeV 2 chosen to give the correct value for f n . This form for 
the quark propagator is based upon studies of the Dyson- Schwinger equation for S(p) using 
a gluon propagator with an infrared singularity, 



g 2 DM= [5, u - k ^ 



k 2 



8ir 4 D6\k), 



(15) 



and a dressed quark-gluon vertex @|J. Quark confinement is characterised by the non- 
observation of free quark states. A pole in the quark propagator in the timelike region 
would be a sign of a free quark state, just as a free bound state is manifest as a pole in a 
n-point Green's function. Thus a sufficient condition for the lack of free quark production 
thresholds is the absence of timelike poles in the propagator. The model quark propagator 
given above is an entire function (except at timelike-p 2 = oo) and so does not have a 
Lehmann representation. This means it can be interpreted as describing a confined particle 
and it ensures the lack of the unphysical singularities corresponding to free quarks in I a ^p. 

The dipole moments are also calculated using a form for the quark propagator developed 
by Mitchell and Tandy f7j, to investigate p-uj mixing. This propagator is given by 



ag(x) = C m e~ 2x + - (l 
x v 

2x) 



av(x 



-(1 

2x 2 



-2x 



(16) 



(17) 



To fix the parameters A = \JlD and C^, a fit to (qq), f^, r n and the ir-ir scattering lengths 
was done. With \(m u + rrid) = 16 MeV a best fit was obtained for A = 0.889 GeV and 



Cm = 0.581 0. This form has a deficiency which can be seen in it's failure to correctly 
model the behaviour of as away from x = 0, in the massless limit with dynamically broken 
chiral symmetry [|J . The large value for the mass of the quark used is related to the deficiency 



of the propagator mentioned above [I3J 



E. Rho-meson Amplitude 



The dominant Bethe-Salpeter amplitude for the rho meson is given by [14 , 15] 



r.„ (M = ,( 7 ,_Jfe)^M (18 ) 

where k is the relative momentum of the quark and anti-quark, p is the momentum of the 
rho meson and I and \x are flavour and Dirac indices respectively. This form ignores other 
allowable Dirac structure in the vector meson Bethe-Salpeter amplitude and so introduces 
errors of the order of 10% ||. Using the quark propagator defined in equations (12) and (13) 
and a Ball-Chiu quark-photon vertex (see later), Chappell uses the following approximate 
form for T p (k,p) [|] 



with ax = 0.38845, ci2 = 0.01478, a = 2. The values for the parameters were found by fitting 
to the experimental values of f p and g p7T7T ||. Pichowsky and Lee, also using the quark 
propagator defined in equations (12) and (13) and a quark-photon vertex of the Ball-Chiu 
type, use an identical form for T p , given by 



r '*»> = e -"* + iTW (20) 

Fitting to the experimental values for f p and they obtain values for the parameters of 
ay = 0.400, by = 0.008 and cy = 125.0. Using the quark propagator defined in equations 
(15) and (16) Mitchell and Tandy use a form for the amplitude T p given by 

r p (k,p) = e- k2 / a \ (21) 



where a = 0.194 GeV. The normalisation for the rho amplitude is fixed by [|10 



P, - ^) = N c tr D J ^ d -^Y pa {k, P )S{kJ)Y p ,{k, P ) (22) 

+iV c tr D J -0^S(k + )T pa (k,p)^^-T pf3 (k,p), 

where k a = k + 2p. This condition, along with the fact that the quark-photon vertex given 
below obeys the Ward Identity, ensures that S(q 2 = 0) = 1, i.e., that the rho has unit charge 

Using the following transformation properties: S T (—k) = C^S(k)C; T^(—k,p) = 
-C^ a {k,p)C and rjf+ cp T {-p, -q) = -C^ c+CP {p,q)C, where C = 7274 is the charge 
conjugation operator, and the Ward-Takahashi Identity, one can show that current conser- 
vation holds ||], i.e., 

q„Ia„p = 0. (23) 

Given a description for the quark-photon vertex, the dipole moments of the rho can now 
be calculated using equations (12), (13) and (18); (12), (13) and (19), and (15), (16) and 
(20). 

F. Quark-Photon Vertex 

The quark-photon vertex also satisfies its own Dyson-Schwinger equation, but solving 
this integral equation is difficult. Despite this, a realistic ansatz for the vertex function has 



been developed by constraining its form using certain criteria |Tj,|T8j. The quark-photon 
vertex ansatz thus obtained is given by 

rf +CP (p,g) = rf(p,g) + rf( P ,g). (24) 

The Ball-Chiu vertex, r PC has the form 0] 

rf(,^^±^7, + 

- z(B(p 2 ) - B{q 2 )) . (25) 



This vertex ansatz is completely described by the dressed quark propagator and satisfies 



both the Ward-Takahashi and Ward Identities, is free of kinematic singularities as q 2 



P , 



transforms correctly under appropriate transformations and reduces to the perturbative 
limit. To ensure multiplicative renormalizability Curtis and Pennington added a transverse 
piece to the Ball-Chiu vertex This term has the form given below 



with 



(26) 



dip, q) 



1 

p2 _|_ g2 



(p*-qZ) 2 + (M 2 (p 2 ) + M 2 (q' 



(27) 



M{p 2 ) 



B{p 2 ) 
A(p 2 ) 



G. The Magnetic Dipole Moment 

The dipole moment can now be calculated. The colour, flavour and Dirac traces are 
performed for the integral J aft a, given in appendix A, and the charge and dipole terms 
projected out using the projection operators outlined above. The integrals are performed 
and the dipole calculated using 

fi = lim M(q 2 ) 
<? 2 — 

= (28) 
The calculation is performed in the Breit Frame where 

g = (0,0,Q,0), 

p= (0,0,3/2,^ + Q 2 /4), 
p'= (0,0,-3/2,^ + 32/4). 

A transformation to 4-dimensional hyper-spherical coordinates is carried out so that 
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k = \k\ (sin/3 sin 6 cos </>, sin /3 sin 6*sin0, sin/3 cos 9, cos/?), (29) 

with 

/fOC fTT f"K r2iv 

d 4 k= k 3 dk sin 2 pdp smOffl d(j). (30) 
Jo Jo Jo Jo 



There is no dependence in the integrand so the integral is performed analytically, con- 
tributing a factor of 2n. 

The integral is then performed numerically using Gaussian quadrature methods to obtain 
the results given in Table |. The results obtained using the framework of Chappell, and 
Pichowsky and Lee, are within 10% of each other which is indicative of the similarities of 
the approaches. The perturbative result from ]TT|] , and the non-relativistic SU(6) result are 



given for comparison. The later result is simply f22 |, 



^ = 2< (31) 

In equation (31) the quark mass is to be interpreted a constituent quark mass of say 340 
MeV [pB| , which gives fi p + = 2.26 e/2m p . Thus one could say that the effect of QCD as 
interpreted in our calculation is to convert m q from the current quark value to an appropriate 
constituent quark value. Note that the improper use of a current quark mass in equation 
(31) gives the large value fj, p + = 152 e/2m p . 



H. The Electric Dipole Moment 

CP violation is a feature of the weak interaction which is not well understood at present. 
It will, in principle, give rise to electric dipole moments for quarks and hadrons. In calculat- 
ing hadronic electric dipole moments, the question of the effects of confinement and QCD 
in going from the quark EDM to the hadron EDM have not been studied. Here we study 
these effects for the rho meson. 

The magnetic dipole of the rho meson comes from the term in the Salam-Delbourgo 
vertex function 
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(9a»qp -gpy.q a )M. 
This results in a term in the pp'-f matrix element that looks like 

e(p).F.e(p')M, 

where 

F a /3 = dpAa — d a A/3 

is the electromagnetic field tensor. To find the vertex function term that contributes to the 
electric dipole moment the familiar correspondence is used: 

M(q 2 )^V(q 2 ), (32) 
where T>(0) = d, the electric dipole moment in units of e/2m. Thus as 

the term in the pp7 vertex that defines the rho EDM is 

e a ^q v V(q 2 ). (33) 

Once again the electric and dipole terms will need to be projected out. The vertex 
function is given by 

= d a "(p) [P,(-g ua E(q 2 ) + q u qMq 2 )) 

+(9»vq* ~ 9^q v )M{q 2 ) + e vaw> q? V(q 2 )] df{p'). (34) 
Using a projection operator of the form 

= - (t^) d W (p)e^\d^'(p'), (35) 



the result 



P^V^ Um V° V(0) (36) 



follows. 
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I. Quark Contribution to the Electric Dipole Moment 



To introduce an electric dipole moment in the impulse approximation the qq*y vertex 
needs to be modified to introduce a CP-violating term. Thus the following vertex is used: 



r>, q) = rf +CP (p, q) - iD qlb a^{p - q)", (37) 

where D q is the EDM of the quark. This form for the qq'y vertex still ensures that both 
the Ward-Takahashi and Ward Identities are preserved, is free of kinematic singularities 
as q 2 — > p 2 and reduces to the perturbative limit. It also still preserves multiplicative 
renormalizability. However we have ignored any structure of the CP-violating term which 
may be generated by strong or electromagnetic interactions. 

Following the same procedure as outlined above, the impulse approximation to the EDM 
is given by 

M fi = s a (p)e^q)e^p')ei:% (38) 

where the integral J^™ is given by 

= (- X ) / ^tr CFD [r%(k;p-q)S(k- + )zr,(k ++ ,k- + ) 

xS(k ++ )I*(k + q/2;p)S(k + -)]. (39) 

The expanded form of the integral is given in appendix B. The EDM can now be 
calculated. Once again the Dirac traces are performed and the charge and electric dipole 
terms projected out using the projection operators given above. The integrals are performed 
and the dipole calculated using 

d = lim V(q 2 ) 
= limP^/J. (40) 

The integral is again evaluated in the Breit frame. 

The dipole interaction between the quark and the photon is governed by the term 
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Duii^a^q^uA^ + D d dl^ u q v dA^ (41) 

where the D u and Dd are the electric dipole moments for the u and d-quarks respectively. 
In many models D oc Q for quarks, so the following form is used, 

D{ u ,d} = Q{u,d} d q . (42) 

The perturbative result has been calculated and it gives: 

rfP ert - = 3m q d q in units of e/2m. (43) 

The quark electric dipole moment has been left as a free parameter and so an estimate 
of the size of d q is needed. The quark electric dipole moment can be generated in various 
models. One of these is the Weinberg model of three Higgs doublets |[24|| . In this model 
the quark EDM comes about due to neutral Higgs exchange with the quark dipole moment 
given by |25| 

^ = ^™,E^/4-'-#. (44) 

where the mm are the masses of the neutral Higgs and X- and Y( are mixings of the neutral 
Higgs. This gives a quark dipole moment of the order of 10 -24 e.cm. Thus we obtain the 
results for the rho electric dipole moment, in units of d q , listed in Table |J. Again the 
perturbative result is included for comparison. For completeness, it is noted that the non- 
relativistic result is given by d p = d q e.cm. The results from the Chappell, and Pichowsky 
and Lee, amplitudes are again seen to be very close, as was the case for the magnetic moment. 



III. THE BAG MODEL 

There are purely phenomenological models that have been used over the years to de- 
scribe the properties of hadrons. One successful and popular model which incorporates 
many of the features of QCD is the MIT bag model fl2"E| , |2"T| l . The MIT bag model includes 
relativistic, point-like, confined quarks and provides a phenomenological description of the 
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non-perturbative gluon interactions. The simplest form of the MIT model will be used, with 
one-gluon exchange and centre of mass motion ignored. While these results, at least for 
the magnetic moments of baryons, are well known we discuss them here to provide another 
point of comparison for the results of section [TJ. 

The starting point for the bag model is the assumption of free quarks confined within 
a spherical volume of radius R. Solving the Dirac equation inside the sphere results in a 
wavefunction given by [p7,2q 

( I w+m „■ ( xr\ A/ \ 



N(x) 



An 



\ 



-1(7 



jo (i) X 
• rji (f ) X 



(45) 



where w(m, R) = y m 2 + x 2 /R 2 , j and j\ are spherical Bessel functions and x is the quark 
momentum in units of 1/R. The normalisation, N(x), is given by 



N-\x) = R*f (x) 2 ™ iw - 1/R)+m « /R 



w(w — m) 



(46) 



and the eigenfrequency x by 

tanx 



x 



1 — m q R — (x 2 



(m q R) 2 ) 1/2 '' 



(47) 



arising from the boundary condition: 



— 27 ■ rip = ip, at r = R. 



(48) 



To confine the quarks a "pressure" B is introduced which ensures energy-momentum con- 
servation and confinement. Thus the total energy of the bag becomes 

E{R) = J2iNi (m% + x 2 /R 2 ) V ' + B^irR 3 , (49) 

where Ni is the number of quarks or type i inside the bag. The stability of the bag also 
implies that 

dE(R) 



dR 



0. 



(50) 



All calculations will be carried out in the rest frame of the bag and so E(R) will be equated 
with m p . Equations (|47D , ( ^9|) and (|50|) will then be used to fix the variables x, R and B. 
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A. Magnetic Moment in the Bag Model 



The magnetic moment is given by 



1 f 



where the summation is over the quarks. This gives [pq ,[23 



= HqEi^iQi, (52) 

where fi q , the magnetic dipole moment of the quark, is given by 

e 1 4wfl + 2m g fl - 3 

^ ~ 2m, 3 m " 2( Wj R) 2 - 2^ + m q R' [b6) 

Using the standard spin-flavour wavefunction of the rho yields the results, listed in Table [TIT], 
and discussed in the final section of this paper. 



B. Electric Dipole Moment in the Bag Model 

Using the single particle wavefunction from above, and the P and CP violating electro- 
magnetic current 

J^u = -iDg^a^ip, (54) 
where a^ u = ~ [7^, j v ] it can be shown that 

J^F^ = -Dqip&ip ■ E + terms involving B, (55) 
where F^ u is the standard electromagnetic tensor. Therefore |29 
d P = f (P T I E D q j><riJ> \p T) 

Jba S quarks 

^r^^n'O^i s^h). (5 6 ) 



quarks 
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Thus, again using the rho spin-flavour wavefunction, the following formula for the EDM of 
the rho is obtained, 

d p = N 2 jf * dr r 2 (j* + ij\ 2 (^^j , in units of d q e cm. (57) 

This leads to the results given in Table [TV] and discussed below. 

IV. DISCUSSION 

Using a model form for the quark propagator obtained from the Dyson- Schwinger equa- 
tions, two different phenomenological rho-meson amplitudes fitted to f p and g p7Tn and one 
used to describe p-uj mixing, and a dressed quark-photon vertex, the rho dipole moments 
have been calculated. 

All three non-perturbative models used, agree qualitatively, with the amplitude and 
propagator defined in |J yielding the largest value for fi p and that of [0] the smallest. The 
Bethe-Salpeter amplitude of QTUJ resulted in the largest value for d p , with that of |7| yielding 
the smallest. Thus the values calculated for the rho magnetic and electric dipole moments 
are fairly robust when it comes to slight changes in the forms for both the quark propagator 
and the rho Bethe-Salpeter amplitude. 

As seen above the inclusion of quark confinement gives greatly enhanced values over 
those obtained using perturbation theory. The non-perturbative results for the magnetic 
dipole moment are ~29%-51% bigger than the perturbative result. While the electric dipole 
moment showed an even greater increase, being ~63-78 times larger than the value obtained 
with perturbation theory. d p is linearly dependent upon the quark dipole moment and so a 
different model for the quark moment could change the value for d p substantially. Conversely, 
an experimental measure for d p would put a limit on d q and possibly rule out some models. 

The increase in the dipole moments is due to the inclusion of quark propagators that 
ensure confinement and a more accurate modelling of the momentum distribution and finite 
size of the rho meson. The kernel of the Bethe-Salpeter equation involves the dressed quark 
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and gluon propagators and so includes the effects of the quark-gluon sea. The amplitudes 
that have been used here can thus be said to model these effects phenomenologically. 

These calculations extend the process of using the Dyson-Schwinger/Bethe-Salpeter 
framework to obtain hadron properties, to the rho meson. 

The rho dipole moments have also been calculated in the MIT bag model. This phe- 
nomenological model also includes confinement via the introduction of the "pressure" term 
B. The bag is also assumed to incorporate the effects of the non-perturbative gluon ex- 
change and interaction. This relativistic, confining model also yields values for the dipole 
moments greater than that of perturbation theory, as was the case for the results using the 
DSE/BSE framework, i.e., the results from both approaches agree qualitatively. The mag- 
netic dipole moment is about ~10% bigger than the perturbative result, smaller than any of 
the results obtained in the impulse approximation. The electric dipole calculation was ~6% 
bigger than the largest result using the Dyson-Schwinger/Bethe-Salpeter equations and ~83 
times bigger than the perturbative result. Again the inclusion of confinement has led to 
an enhancement in the dipole moments of the rho. Indeed, one could suggest that for the 
calculation of the magnetic and electric dipole moments, the bag model gives a reasonable 
estimate of the effects of confinement. 

The rho meson has been used in this study, not because there is any realistic possibility 
of measuring the electromagnetic moments we have calculated, but to provide some insight 
into the effects of QCD and confinement on calculations of these moments in the hope of 
obtaining insight into the effects these aspects of quark dynamics may have on estimates of 
the electromagnetic moments of hadrons • 
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APPENDIX A: INTEGRAL FOR THE MAGNETIC DIPOLE MOMENT 



With the definitions given above the integral J aw g becomes (after colour and flavour 
traces) 



2lN r 



d A k 



T p {k)T p {k + q/2) tr D 



7/3 #-+7^ ft++la h- T i + 7/3 #-+7„7 a T 2 



+7/3 /^-+7a ft+-k 0+f ,T 3 + 7^ J6_ + j£ ++ -y a k 0+IJi T 4 
+7/3 ^_+7a h-ko+^Ty + 7/3 ji ++ -f a k 0+lx T 8 

+ll3l^la ft+-Tg + 7,37m ^++7«^10 

+7/3 ^++7a ft+-k 0+fJ ,T n + 7/37^^0+^712 

+7/3 ^-+7« ft+~k 0+II T l3 + 7/3 ^ ++ 7 a £:o+/ t Ti4 

+ 7/3 /^++7a ^+-^0+^^15 + 7/37^0+^16 



(Al) 



where the terms with an odd number of 7 matrices have been dropped. The 7\_ 16 are 
defined as follows: 



7i = i 



i(K + V A )a v (k 2 _ + )a v (k: 



T 2 = -i{V 1 + V 4 )a v (k 2 _ + )< 
V 5 )a v (k 2 _ 

MV2 + V> y (A; 2 > 5 (A; 2 + K(A; 2 _), 

T 5 = 



T 3 = 2^ 2 + (\/ 2 + 



)<7v(fc+_), 
) C T 5 (4 + ) C 7 5 (4_), 

)a y (A; 2 



V+ 



loy(fc + _), 



-, = 2zfci + (y 2 - V 5 )(Jy(/c 2 . + )(Jy(/c 2 r+ )(Ty(A; 2 ._), 

T 6 = -rf + (V2 - V> y (A; 2 + K(A; 2 + K(A; 2 _), 

T 7 = -2^^-+)^++)^+-), 
T 8 = -2V 3 a v (k 2 _ + )a v (k 2 ++ )a s (kl_) : 
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T 9 = 


/it- , t r \ /j 2 \ /j 2 \ /j 2 \ 

-i(Vi + V 4 )a s (ki + )a s (k 2 ++ )a v (k 2 + ^), 


T w = 


-i(Vi + V 4 )a s (kt + )a v (k 2 ++ )a s (k 2 + _), 


Tn = 


-2z(T/ 2 + \/ 5 )a 5 (A;i + ) ( T5(A;; + ) f 7y(A;;_), 


T 12 = 


-2ik 2 ++ {V 2 + V 5 )(T5(A;l + )cry(A;; + )cT 5 (A;;_), 


Tl 3 = 


-2i(V 2 - V 5 )a s (k 2 _ + )a s (k 2 ++ )a v (kl_), 


T u = 


-2i{V 2 - V 5 )a s (k 2 _ + )a v (k 2 ++ )a s (k 2 + _), 


Tl 5 = 


-2V 3 a s (k 2 _ + )a v (k 2 ++ )a v (k 2 + _) : 



T 16 = 2V 3 a s (k 2 _ + )as(k 2 ++ )a s (k 2 + _), 



with, 

Vl ~ 2 ' 

- ! (g(g ± ) - g(g + )) 
k 2 — k 2 ' 

2rf(A; ++ ,A;_ + ) 
y -(A(^ + )-A(A; 2 + )) 
5 2d(k ++ , k_ + ) ■ 



APPENDIX B: INTEGRAL FOR THE ELECTRIC DIPOLE MOMENT 

The integral J^™ becomes (after colour and flavour traces) 

7/3 ^-+7 M ^++7a + 7/3 ^-+7 M 7aT 2 

+7/3 ^-+7<* j£ + -k 0+II T 3 + 7/3 ^ ++ 7 Q A;o+ M T4 
+7/3 ^-+7a j£+-k +»T 5 + 'y^ako+^Ta 

+7/3 ^++7a^0+M T 7 + 7/3 /^-+7a ft+-k 0+fl T 8 

+7/3 ^_+7 5 o- w g p ft++l a T 9 + 7/3 J6-+'y 5 <T lip q p 'y a ft+-T w 
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+7/37m h+laTn + 7/37 M 7a #+-7\ 2 

+7/37^0+^13 + 7/3 jt+^a ft + -k 0+II T 14: 

+7/3 H-+ ^++lak +^,T l5 + 7/3 ^_+7 a ^ + ^k 0+fl T w 

+7/3 ^++7a ft + -k 0+II Ti 7 + 7/37a^0+M T 18 



+7/3750" /ipg P ^++7a ^+-^19 + 7/375^^7^ 
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(Bl) 



Again dropping terms with an odd number of 7 matrices. The Ti_ 20 are defined as follows: 

T x = ifr + V 4 )a v (k 2 _ + )a v (k 2 ++ )a v (k 2 + _) 

T 2 = -ziV, + V 4 )a v (k 2 _ + )as(k 2 ++ )a s (k 2 + _) 

T 3 = 2^ + (\/ 2 + ^ 5 )M*- JM*++K(*+ _) 

T 4 = -2i{V 2 + V 5 )M*- + )a 5 (^ + ) (75 (A;2_) 

T 5 = 2ik 2 _ + (V 2 - V b )a v {k 2 _ + )a v {k 2 ++ )(j v {k 2 + _) 

T 6 = -2iki + (V 2 - V 5 )a v (k 2 _ + )a s (k 2 ++ )a s (k 2 + _) 

T 7 = -2V 3 a v (k 2 _ + )a v (k 2 ++ )a s (k 2 + _) 

T 8 = -2V^ v {k 2 _ + )a s {kl + )a v {k 2 + _) 

T 9 = -V 6 a v (k 2 _ + )a v (k 2 ++ )a s (k 2 + _) 
T 10 = -V 6 a v (k 2 _ + )a s (k 2 ++ )a v (k 2 + _) 
T u = ~i{V x + V 4 )a s (k 2 _ + )a v (k 2 ++ )a s (k 2 + _) 
T 12 = -i(Vx + V 4 )a s (k 2 _ + )a s (k 2 ++ )a v (k 2 + _) 
T 13 = -2ik 2 ++ (V 2 + V 5 )as(k 2 _ + )a v (k 2 ++ )a s (k 2 + _) 
T 14 = -2%{V 2 + V 5 )a s (k 2 _ + )a s (k 2 ++ )a v (k 2 + _) 
T 15 = -2z(F 2 - V 5 )as(k 2 _ + )a v (k 2 ++ )a s (k 2 + _) 
T 16 = -2z(F 2 - V 5 )as(k 2 _ + )a s (k 2 ++ )a v (k 2 + _) 
T 17 = -2\/ 3( r 5 (A; 2 + ) ( T y (4 + ) ( T y (4_) 
T 18 = 2V 3 a s (k 2 _ + )a s (k 2 ++ )a s (kl _) 
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T 19 = -V 6 a s (k 2 _ + )a v (kl + )a v (k 2 + _) 

T 20 = V 6 as(k 2 _ + )a s (kl + )as(k 2 + _), (B2) 
where V\. . . V5 are defined as before and 

Vq = —id g , d q the quark dipole moment. (B3) 
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TABLE I. Magnetic Dipole Moment fi in units of e/2m 



Rho Amplitude/Quark Propagator 




Chappell @ 


3.01 


Pichowsky and Lee 


2.72 


Mitchell and Tandy j§ 


2.57 


Perturbative Result jllfl 


2.00 


Non-Relativistic SU(6) [p^|, constituent quark mass from 


2.26 
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TABLE II. Electric Dipole Moment d p in units of d q e.cm. 



Rho Amplitude/Quark Propagator 


dp 


Chappell [| 


0.743 


Pichowsky and Lee 


0.779 


Mitchell and Tandy 


0.627 


Perturbative Result jn], m q from PJlO[ 


0.010 
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TABLE III. Magnetic Dipole Moment fi in the Bag Model in units of e/2m 



Mass of Quark 




m q from [p|JlO[| 


2.21 


Chappell @ 


3.01 


Pichowsky and Lee [lOf 


2.72 


Mitchell and Tandy @ 


2.57 
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TABLE IV. Electric Dipole Moment in the Bag Model in units of d q e cm. 



Mass of Quark 


dp 


m q from PJlO[| 


0.828 


Chappell @ 


0.743 


Pichowsky and Lee 


0.779 


Mitchell and Tandy 


0.627 
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FIG. 1. The dressed p/ry Vertex 
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FIG. 2. The quark contribution to the dipole moments of the rho, with S the dressed 
quark propagators, the dressed quark-photon vertex and T p the rho-quark vertex. 
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